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When a charge is held static near a constant density spherical star, it experiences a self-force F ae if , 
which is significantly different from the force it would experience when placed near a black hole 
of the same mass. In this paper, an expression for the self-force (as measured by a locally inertial 
observer) is given for an insulating Schwarzschild star, and the result is explicitly computed for the 
extreme density case, which has a singularity at its center. The force is found to be repulsive. A 
similar calculation of the self-force is also performed for a conducting star. This calculation is valid 
for any static, spherically conducting star, since the result is independent of the interior metric. 
When the charge is placed very close to the conducting star, the force is found to be attractive but 
when the charge is placed beyond a certain distance (2.95M for a conducting star of radius 2.25M), 
\ the force is found to be repulsive. When the charge is placed very far from the star (be it conducting 

or insulating), the charge experiences the same repulsive force it would experience when placed in 
' the spacetime of a black hole with the same mass as the star. 
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I. INTRODUCTION 



The electrostatic potential and the fields produced by a static electric charge in the vicinity of a Schwarzschild 
black hole have been discussed in detail in several papers [H, 0, S, 0| ■ These all suppose that an external force holds 
an electric charge e at rest at a coordinate position R — b in the Schwarzschild geometry. This external force should 
balance the electrostatic force on the charge in addition to the usual attractive gravitational force (on the mass 
carrying the charge). When there is no other charge anywhere, the entire electrostatic force on the charge may be 
5-H , regarded as a self-force. 

An explicit expression for the electrostatic potential of such a point charge is given as a summation over multipole 
moments in [l[ . The horizon of the black hole is found to be an equipotential surface and Q shows the electric lines 
£NJ ■ of force everywhere outside the horizon. A closed form expression for a specific potential, V c , was given earlier by 
Copson using work of Hadamard [f|, but it corresponds to a solution with non-zero charge on the black hole. For 
the potential corresponding to physically relevant boundary conditions (namely, zero net charge on the black hole), 
LinetQ observed that a spherically symmetric homogeneous solution had to be added to Copson's potential, yielding 
! a result in accord with 

By employing a strategy similar to Dirac'sQ, of imposing conservation of the stress energy tensor inside a world- 
tube surrounding the particle and then limiting the world tube to the world line of the particle, Smith and WilL7] 
calculated the external force needed to hold this charge at rest. After subtracting the requisite gravitational force, the 
electrostatic self-force is found to be F®^ — e 2 Al/b 3 . 1 It is implicit that this force is always calculated with respect 
to a locally inertial observer at the position of the charge. The potential of [l[ or [3| is essential for understanding 
^ ' this result: with the potential, V c , given by CopsonQ, no self-force would have been found. 

A static electric charge is a singular point source and the potentials discussed so far all satisfy Maxwell's equations 
with the point particle as source. This is quite distinct from the behavior of the so-called "direct" and "tail" fields 
Q, for which the source in Maxwell's equations is not well identified. By contrast, the radiation potential of DiracQ 
satisfies Maxwell's equations with no source. Other familiar potentials which satisfy Maxwell's equations with the 
point particle as source include the retarded and advanced potentials, and their symmetric sum. Ironically, Copson's 
potential in 0] is none of these familiar constructs, nor is the potential of [H, H[. 

Detweiler and Whiting Q demonstrate that, in general, a potential V due to a point source can be written as sum 
of two parts, V = V s + V~ R , where V s is a singular potential (divergent at the particle) which exerts no self force on 
the particle and V R is the regular potential which is entirely responsible for the self-force. In particular, locally near 
the source, V s satisfies Maxwell's equations with the point particle as source, while V R is an homogeneous solution, 
without source. Within a normal neighborhood of the point charge, the singular potential can be obtained from a 
variant of Hadamard's form of the Green's function^. We will need to understand its relation to the potential, V c , 
of Copson, that is, locally. 



1 We use units in which G = c = 1 throughout. 
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A. The Singular Potential 

It can be understood that Copson's potential for this static problem is exactly V because of the relation between 
the construction of V s in Q and Copson's construction of V c using Hadamard's work. In constructing the potential 
we refer to as V c , Copson used the unique, locally- represented, least-singular (hence) elementary solution for the static 
problem, as given in a general formalism due to HadamardQ. By construction, the singularity in V c at the source is of 
as low an order as is possible. In the not necessarily static case, the singular potential V s is, in a normal neighborhood 
of the point charge, the unique locally-determined, similarly least-singular solution of generalized Hadamard formQ, 
which has support on and outside the light-cone (it exerts no force). In this regard, V s is unique, due to its lack of 
(local) support inside the light cones. In the static case, when the light-cone collapses and time derivatives disappear 
in the equation for the potential, the domain of support for V becomes precisely the domain in which V is defined. 
Then, V and V solve the same problem under the same conditions, and hence coincide. 

The result that Copson's potential V c is precisely V s , is indeed compatible with the fact that it yields no self- force. 
Consequently, the correction to Copson's potential given in Q automatically corresponds to V R , the regular part of 
the potential, which provides the self force. 



B. Overview 



In this paper, we consider replacing the black hole by a constant density star. The fundamental fact we will use is 
that, locally, the potential constructed by Copson V c is exactly the singular piece V s , as long as the metric in the 
neighbourhood of the charge is the Schwarzschild metric. Any additional part of the potential at the position of the 
charge will contribute to the self force. Our aim is to calculate the self force of a charge e placed in the vicinity of a 
spherical star in two different situations — depending upon whether the star is totally conducting or totally insulating. 
Since the metric outside the star is the Schwarzschild metric, Copson's potential is indeed V s . To calculate the self 
force, we follow Linetjl] by first finding the homogenous solution of Maxwell's equations V, which should be added 
to V^ s (given by Copson) to give the whole potential. 

In section [ill wc briefly discuss the Schwarzschild black hole metric, and the metric of a constant density 
(Schwarzschild) star, which we show is conformally related to a much simpler metric. In section IIII1 wc obtain 
the solutions to Maxwell's equations inside and outside the star. We start with the singular potential V s outside 
the star and continue the solution inside the star's surface. This reveals that the surface of the star carries a charge 
distribution, a result which is not physically acceptable. We correct this in section 4, where we impose two different, 
physically interesting, boundary conditions to calculate the homogeneous solution V, which in turn is responsible for 
the self force. In section HV A[ we consider an insulating star that does not get electrically polarized. An additional, 
regular field, arising both inside and outside the star, is necessary to cancel the surface charge induced by V s . It is 
this additional field which is responsible for the self force on the particle. We find that the self force F se \f , is repulsive 
and is significantly different from F®^ when the charge is placed very close to the star. In section HVBl we consider a 
conducting star which has no electric field in the interior, but nevertheless requires an additional field in the exterior 
to remove the net surface charge which would otherwise be induced on the star. Here, we find that F sc \f is attractive 
when the charge is placed close to the star (b < 2.95M) and is repulsive when the charge is placed farther than 2.95M. 
We should interpret the attractive force as a consequence of electrical polarization of the star and the repulsive force 
as due to interaction of the fields with space-time curvature. 



C. Discussion 



For calcula ting the self force of a point source in curved spacetime, a variety of regularization techniques have been 
used pj. [qL Hoi lllj. In a case somewhat related to ours, with a point source placed near a thin shell, Burko et al.[l2j 
have explicitly shown how the regularization techniques can be used to calculate the self force. In this paper, we do 
not perform a regularization to calculate the self force, but we very much rely on the regularization details in @, Q- 

We have focused on the extreme density case because of its intrinsic interest. In principle, our approach here can be 
used to calculate the self-force on a static electric charge in the vicinity of any star with a spherically symmetric metric, 
not just a star with constant density. Generally, the resulting expression for the self-force would be very complicated 
and would require numerical evaluation, which is actually what we resort to in the final step of evaluation of the 
self-force in this paper. It would be interesting to investigate whether such results could have any detectable effect 
on neutron star physics and observations. 

Another question that might come to mind is whether the self- force calculated for the static charge F Be \{, has any 
connection with the self force experienced by a freely falling charge on which no external force acts. DcWitt and 
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DcWitt [l3| have shown that the self force on a moving charge consists of two parts, a radiative part, which removes 
the energy from the particle and dumps it into the fields and a conservative part, which pushes the particle off a 
geodesic without removing energy from it. For a slow moving charge, they explicitly show that the conservative 
part of the self-force is independent of the particle's velocity. In the context of the charge placed near a black hole, 
Wiseman|l4[ clearly points out that, the self force F®^ experienced by the static charge (as calculated by Will and 
Smith) exactly corresponds to the conservative part of the self-force experienced by a freely falling charge. It is 
straightforward to see that we can adopt the same interpretation for the self-force calculated in this paper. 



II. THE METRIC 



Consider a constant density star of total mass M and radius R s . The metric outside the surface of the star is the 
Schwarzschild black hole metric (see, for example, [HI): 

ds 2 = -(1 - 2M/R) dT 2 + d fl I/R) + R 2 dfl 2 , R>R S , (1) 



and the metric inside the (Schwarzschild) star (see also [Ta|) is 
ds 2 = 



a- iv 7 ! -aR 2 



dR 2 



dT +m ^i+ Rdn > R<R s , (2) 

[1 — aR z \ 



where a = | y/1 — 2Mj R s and a = 2M/Rg. The quantity M/R s is bounded above by 4/9. Hence the range of a is 
also restricted: 

0<M/i? s <4/9 =S> 3/2>a>l/2. (3) 

For the extreme density case, when M/R s = 4/9, we have a =1/2 and the center, R = 0, develops a singularity (the 
component, Grr, of the Einstein tensor blows up). 

The Weyl tensor for the metric in ([2]) evaluates to zero, which means that the metric is conformally flat. This 
interior metric can be written in the following simpler form: 

ds 2 = s [-([3 2 + ar 2 ) 2 dt 2 + dr 2 + rW] , (4) 

(1 + ar 2 ) 

where 1 = (2a — 1)/ (2a + 1) and the new coordinates are defined as: 



t = T(2a+l)/4, and r = R/(l + y/l - aR 2 ). (5) 

The coordinate r is a monotonic function of R, and r ranges from to some maximum value, r s at the surface of the 
star. We have the following expressions for r s and /3 2 : 

R ° P 2 = \~ 2ar l (6) 



1 + 2M 1 - ar- 

For each constant density star, f3 2 is a constant, and it can take values which range between and 1/2. For the 
extreme density case, we have (3 = 0, and r s = 3i? s /4 = 27M/16, in which case the interior metric takes the following 
form: 

ds 2 = ? \~a 2 r 4 dt 2 + dr 2 + r 2 dQ 2 ] . (7) 

(l + ar 2 ) 2 L J 

There exist coordinate transformations which make (H]) manifestly conformally flat[l(|. Except possibly for (3 = 0, 
those coordinate transformations mix the spatial and temporal coordinates. Our focus is restricted to electrostatics 
in this paper. Since Maxwell's equations are conformally invariant, it will be sufficient to work with the conformal 
metric evident in ((4]), obtained by ignoring the conformal factor 4/ (l + ar 2 ) . Thus we will use the metric: 

ds 2 = [-(ft 2 + ar 2 ) 2 dt 2 + dr 2 + r 2 dft 2 ] , (8) 
when solving for the electrostatic potential in the following section. 
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III. ELECTROSTATICS 



Maxwell's equations, which govern the theory of electrodynamics, can be written as 

V^F"" = iirJ", (9) 

where = d^A v — d^A^. In terms of the coordinate derivatives, they can be presented in the following explicit 
form: 

dnW=9<f p 9 w Fp*) = 47rV=ff J". (10) 

It is evident that Maxwell's equations are conformally invariant: remains invariant under a conformal transforma- 
tion of the metric provided the source term in (|10|) is defined so that it, too, is invariant under such a transformation. 2 
For electrostatics in a static geometry, Maxwell's equations in the Lorentz gauge (V M ^4 M = 0) can be reduced to a 
single equation: 

il: x —^<r<),.\„> = ^4~g J°. (11) 

In this paper, we consider a charge, e, on the z-axis at the coordinate position R = 6, located outside a static, 
spherical star. The resulting potential Aq, which henceforth we denote as V(r,9), will be axially symmetric about 
the z-axis. 



A. Outside the star 



To find the potential outside the star, we will solve Maxwell's equations, which have been reduced to the single 
equation in (jlip . In a static, spherically symmetric geometry, this takes the form: 



d r (V=99 rr 9 00 d r V) + de(V=gg ee g 00 d e V) = iir^J , 

where 

27iV^J° = e6(R - b)6(6) = eS(R - b) sin 66(1 - cos6>) 
= e6(R-b) sm6J2t [Q + §)P/(cos0)] . 

Using the metric in (fTJ), equation (fT2]) can be written as 



suit 



d f^dV 



1 



1 g ( . 9 dV 
dR\" dRj ' (1-2M/R)sm6d6\ m 06 



2e6(R-b) sin(9^] (£+-)P £ (cost 



(12) 



(13) 



(14) 



We can decompose the solution in terms of Legcndrc functions as 

V(R,6) =^2[Vt(R)P t (coa0)] 



(15) 



and hence simplify (fT4"]) to get 



OR 



OVA e{£ + l)V e 
OR J + (1 - 2M/R) 



= e(2i + l)6(R-b). 



(16) 



Copson[3f constructed a particular solution, V (R,6), to these equations directly from Hadamard's elementary 
solution[5j . Even though the solution he constructed was in the context of the charge in the presence of a Schwarzschild 



2 With remaining invariant under a conformal transformation, A v is also taken to be invariant. Then, in terms of A v , the invariance 
of the other of Maxwell's equations, V[ a .F^„] = 0, under a conformal transformation, follows trivially from the identity, R^^^j = 0, for 
the Riemann tensor, which holds in any metric. 
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black hole, his solution will still be a particular solution for the case where the black hole is replaced by a spherical 
star, because the exterior metric of the spherical star is the same as the black hole metric. Copson gives 

V°(R,0) = i - (b~M)(R-M)-M 2 cos9 ^ 

bR ^J(R - M) 2 - 2(R - M)(b - M) cos 9 + (b - M f - M 2 sin 2 9 

A calculation similar to that in Q explicitly shows that this potential does not contribute to the electrostatic self 
force. It has been known for some time [9( that the field produced by any point source can be decomposed into 
two pieces, a locally-determined, singular field which does not affect the particle's state of motion and a regular field 
(which is a solution to the homogeneous field equations in the neighbourhood of the particle) that does; that is, 
V = V s + V R . The singular potential V s has nothing to do with the self-force. It is the regular potential V K which 
is entirely responsible for the self- force. It was explained in the introduction that V s (R, 9) = V G (R 1 9). 

The singular field V s is an entirely local construct. The boundary conditions and external sources play a part only 
in determining the regular field V R . The regular field interacts with the source and makes it experience a force. It is 
crucial in determining the path the charge actually follows. The external field (due to other sources) is in general a 
part of the regular field. Thus, the force experienced by the point charge is a sum of the force due to external field 
and the self-force. In the absence of any external sources, the force due to the regular field alone corresponds entirely 
to the self-force. 

We can decompose the singular potential (fT7|) in terms of Legendre functions as in (jT5J) : V s (R, 9) = 
J2i [Vp(R) Pi (cos 9)]. Using the results of [l|, the multipolc moments corresponding to this potential are given 
by: 

1 f t (b)g t (R) - &fc, for R s < R < b 



where 



g t (R) = { 2,ai tw M ' ( R - 2M )^£ - !) fOT 1 * (20) 
WK > \ 1 for£ = 

are the two linearly independent homogeneous solutions to (fH)]) . Note the £ = Q results: fo(R) = 1/R and go(R) = 1. 
Otherwise, as R^oo, f e (R)^l/R^ +1) and g e (R)^lxR e . 

In the next section, we will calculate the regular field V R that should be added to the singular field V s in the 
vicinity of the charge, so as to satisfy the appropriate boundary conditions. Before that, we solve for the continuation 
of V s inside the spherical star. 



B. Inside the star 



We assume throughout that the spherical star has a constant density as discussed in section ITT] Since Maxwell's 
equations in four dimensions are conformally invariant (as long as the source is conformally invariant), the solution 
inside the star can be obtained by using a conformally related metric and a conformally invariant source in Maxwell's 
equations. Thus, we will use the conformal metric, ([5]), when we solve Maxwell's equations inside the star. The 
interior potential V(r, 9), can be found from Maxwell's equations by using the interior metric ([5]) in (jT2J) : 

• ft\—( 7 ' 2 dV \ 1 1 d t dV\ 

~ Sm [dr \ (/? 2 + ar 2 ) !fr J + (f3 2 + ar 2 ) ^9d9 \ Sm ~d9 J 

We assume that there is no charge inside the surface of the star. Then, J°s/—g is the charge density on the star, 
which can be expressed as cr(9) sm9d(r — r s ), where <j(9) is the surface charge density. 

By axial symmetry, we can decompose the solution as V(r, 9) = ^2 e Ve(r)Pg (cos 9). We can also decompose the 
charge density as a(9) = Y>i<j l Pi(cos9). Now, (J^J) reduces to 

= a l 8(r-r s ). (22) 



= 4ttJ°V^. (21) 



d_ 

dr 



dV, 



1(1 + 1) 



(f3 2 + ar 2 ) dr J (f3 2 + ar 2 ) 



■Vt 
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The homogeneous solutions to (|22|) are obtained in terms of hypergeometric functions F([a, b], [c], For each £, 

there are two linearly independent solutions for V^(r), namely hg(r), kg(r): 

Mr)=< ;-^(l + ^) 2 ^([l' 2 + ^[^l].-^ 2 ) for /3 7^ ^ 

{ar)~ l for /? = 0, 



and 



r 2 "l 2 F([|,l-^],[i-£],-^r 2 ) for /3 ^ 
(ar) f+1 for ,9 = 0. 



Note that /io(r) = 1, is a constant. When (3 ^ 0, we see that, for small values of r, hg(r) ~ r and fc^(r) ~ r — ^ 
Since, for (3 ^ 0, we want 14 to be well behaved at the origin, it should take the form of hi(r). Away from r = 0, 
he(r) is well behaved in the limit (3 — > 0, so we propose that Vg should also take the form of /if (r) for j3 = 0. 3 

To find the precise form of the potential Vi(r) for £ = and /? = 0, we impose the requirement that there is no 
charge at the center of the star r = 0. To examine the charge at the center of the star at r = 0, we consider a small 
spherical Gaussian surface at r = e around the center. When we integrate the field due to this potential over this 
small Gaussian surface, the charge inside it is proportional to (see appendix |A"|) 

r 2 dV (r) 



f3 2 + ar 2 dr 

Only the £ = term would survive because the integral would involve a term of f^(cos#) integrated over the surface. 
Requirement that the charge at the center should vanish implies that for f3 = 0, the potential Vq(t) takes the form of 
ho(r) which is a constant. There is no ambiguity. 

Since we already have the singular field V s outside the surface of the star from ([T7| and (TT5)) , we can now extend it 
to the interior of the star. Continuity of the vector potential A M at the surface of the star along with the coordinate 
transformation that takes the interior star metric from @ to (0} implies that 

dT 

A t r=rs = — A T \ r=Rb . (25) 
at 

From ([5]), we have dT/dt = 4/ (2a + 1). Putting everything together, we sec that, for (3 ^ 0, we can write the singular 
potential V s inside the star as 

VKr) = 1 ^r T A ^ (26) 

(2a +1) [ f (b)g (Rs) - ^5 for t = 0. 

For (3 = 0, we have a = 1/2, which gives dT/dt = 2, in which case 

( fe(b)ge(R s )(^Y for I? 

1 f (b)g (Rs) - #t for£ = 0. 



IV. SELF FORCE OF THE CHARGE 

The static potential we have just constructed has a discontinuous derivative at the surface of the star, implying a 
surface charge density, a, there. In this section, we shall deal with a in two distinct ways, through a potential V , that 
will actually be responsible for the self-force, which we evaluate in the two specific cases we consider: 

• when the star is insulating, cannot be polarized, and has no free charges whatsoever; continuity of the gradient 
of the potential across the surface will ensure no charge is required there, 

• when the star is highly conducting, with zero net charge; free charges arrange themselves on the surface so that 
there is no electric field inside the star. 

In each case, we obtain the regular potential V^, by imposing suitable boundary conditions on the surface of the star. 



3 We thank a referee for clarifying the argument required here. 
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A. Insulating Star 

The singular field V s constructed in section Hill demands the existence of a charge density on the surface of the star. 
Knowing the field outside the star from (TT5)) and inside the star from (|2"rj)) . we can apply Gauss's law (see appendix 
|A"|) to find this charge density. Writing the area element on the Gaussian surface sin 9d9d(f>, as dtt, and following (|A2[) 
from the appendix on Gauss's law, we obtain: 



4ira(9)dn = - R 



2 dV s (R, 



dR 



dfl 



r 2 8V s (r,9) 



R=R S 



/* 2 



dr 



(28) 



For the extreme density star ((3 = 0), we can re-express (|28j) in terms of the multipolc moments of the charge density 
and the potential using equations (fP8|) and ([27]) : 



47TC7 = 



R 2 f e (b)g/(R S 



21 



M 
b 



ft(b)gt(Rs) for^O 

for £ = 0. 



(29) 



Here, <7° corresponds to the net charge on the surface of the star. It is clear that, in general, higher moments of the 
charge distribution are also non vanishing. 

Since an insulating star should not have any charge distribution on its surface, we have to add a regular potential V R 
to the singular potential V , so that the full potential satisfies the boundary condition of zero charge on the surface. 
Clearly, V R is the (otherwise homogeneous) potential obtained by a charge density of —cr(9) on the stellar surface. We 
call this potential V . As usual, we shall decompose this potential into multipole moments: V(R, 9) = T I gVi(R)Pi(cos 9). 

As we already know from the two linearly independent vacuum solutions for Vi (R) outside the star's surface are 
f e (R) and g e (R) from dHJ) and ([20|). For large R, f e (R) -> R-^+V and g e (R) -> R e . Since we require V(R,9) -> 
for large R, we need Ve(R) = Cefe(R). The Ce are constant coefficients whose values are yet to be determined. Before 
finding these coefficients, we must first find the form of the potential Vg(R) inside the star. 

By requiring Vg(r) inside the surface of the star to impose no source at r = 0, we can fix its form everywhere 
inside the star. For the extreme density star ((3 = 0), we require Vi(r) = Dg/r l for all £, where the Dn are constant 
coefficients. Imposing the continuity of Ve across the star's surface, along with the coordinate transformation which 
takes the metric from ([2]) to (j4j), gives a relation between the Ci and Di coefficients: 



BT 

D e = —C e fe(R s )ri, 



dT 1 
and D Q = -^-Cq — . 

at it, 



(30) 



Since the charge on the surface is —a{6), we can use Gauss's law (as shown in the appendix) to write an equation 
analogous to 



4TT<r(9)dn = - R 



2 8V(R,t 



dR 



dn 



dV(r, 9) 



R=R S 



p 2 + CUT 



dfl. 



(31) 



For the extreme density case ((3 = 0), we can re-express pip in terms of the multipole moments. 

f -R 2 C e f/(Rs) - 2g-Ctf e (R.) for i + 
- Ana = < 

[ C for £ = 0. 

Now (|29| and ((32)) can be solved for the coefficients Cg. For £ = 0, we simply have Co = eM/b. For £^0, we find: 



(32) 



e 



2£ 

R 2 Je'(Rs) + — MRs 



2£ 

R 2 s Mb)g/(R s ) + Mb)gt(R a 

ar s 



(33) 



Evaluating the self force 



The regular potential in the vicinity of the charge has been completely determined: Vi(R) = Cefe(R). The self 
force experienced by the charge e is now given simply by Lorentz force law where E is the electric field 
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measured by a locally inertial observer at rest on the charge (i? = b, 8 = 0). Since the electromagnetic stress tensor 
is a rank 2 tensor, under a coordinate transformation, F a/3 = F M „ {dx^ / dx a ){dx v / dx 13 ) , where we take (//, v) to 
correspond to Schwarzschild coordinates and (a, /?) to correspond to the locally inertial coordinates, which we denote 
by {1Z, T, O, <!>}. The diagonal form of metric in Schwarzschild coordinates, implies a simple coordinate transformation 
to the locally inertial coordinates, dx^/dx * = $aQI vWmA)- The electric field in the locally inertial coordinate system 
is then E x = F 0i /y/\g 00 \ gu. 



E-r — - 



dV(R,i 
dR~ 



E B = - 



1 1 dV(R, 6) 
Ryfl- 2M/R 90 ' 



(34) 



On the z-axis (8 = 0), we have Eg = 0, because, for all £, dPe(cosd)/dd\ g=0 — 0. The force experienced by the charge 
is therefore in the radial direction. Since Pt (cos 0)\ s=o = 1, the self force can be summed up over different components 
of £, to give 



F sc \i = ^2 F* el{ = -e C t 



dfi(R) 



dR 



R=b 



Using (|3"3"j) to substitute for Ct, we get 

e 2 M 



Fsdi — —- 



/o'W + e 2 £ 



R 2 f e (b)gt'(R s ) + (2£)ft(b)g e (R s )/(ar s ) 
R 2 J/(R S ) + (2£)fi(R s )/(ar s ) 



h\b)- 



(35) 



(36) 



With a Schwarzschild black hole, as in [7|, in place of the star, the self force experienced by the charge is just 
^seff = e 2 /b 2 [M/b). We now explicitly compute F se if from |36|) and compare it with F^$. From section ITT| we know 
that for the extreme density star, R s = 9M/4, r s = 27M/16 and ar s = 8/27 AI. We can expand the functions ft and 
ge in (|3"6")) to get an explicit expression for F* el{ . Each -F/ olf can be written as 



F , f = — 

self ^2 



~M~ 










(t) 



(37) 



where we have factored out e /b [M/b] from each -Ff elf , to show its relationship with F^. For convenience, we define 
x = M jb, and give the exact form of rf (x) for £=0, 1, 2: 

rf(x) - 1, (38) 



81 



2x 4 (271n3 - 28) 



In [1 - 2x] - x + x' 



- In [1 - 2x] + x + x 2 



(39) 



r] 2 (x) 



945 



2x 6 (5671n3 - 620) 



3x z 
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-x — 3 



In [1 - 2x] 



In [1 - 2x] + .x 3 + 3x 2 - 6x 



6x 2 + 3x 



(40) 



Note that x can take values in the range (0, 4/9), where x = corresponds to the charge being at infinity and x = 4/9 
corresponds to the charge being on the surface of the star. It turns out that, for every £, r] e (x) is positive and for 
large £, r\ (x) — > 0, for all x in the specified range. In Fig[TJ the convergence of r] e (x) is plotted for various values of 
x. When the charge is close to the surface of the star (a; close to 4/9), we must sum over many £'s to get an accurate 
result for the self force. 

As an example, if the charge is placed at b = 3 M (x = 1/3), we can numerically calculate F se \{ to be 



e 
IP 



M 

T 



l. 



31.8 ■ 



M 

T 



(41) 



The self force experienced by the charge at b = 3M, in the presence of this star is 31.8 times stronger than the force it 
would experience if the star were replaced by an equally massive black hole. As x — > (b — > oo), we have F sc if — > F™: 
a charge placed far away from this star, would feel the same force as it would feel when the star is replaced by a black 
hole of equal mass. But as we move the charge close to the star x — > 4/9, -F S oif becomes orders of magnitude greater 

0.43, we have -F se if 



than -F S gff , for example, at x 



36777 F™. 
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B. Conducting Star 

So far, we have considered an insulating star which cannot be polarized. Now, we turn our attention to a conducting 
star which has free charges on it, but no net charge. Hence there could be non zero induced surface charge density. 
We call this surface charge density a (8) (which can be decomposed as a (0) = [5^(cos0)]). This star would be 
polarized in such a way that the field inside it vanishes. The interior metric of the star plays no role in determining 
the field outside the star, only the size of the star R s matters. As before, V = V s + V is the total potential, V s 
is the singular piece which exerts no self force and V (homogenous in the vicinity of the particle) is obtained by 
incorporating appropriate boundary conditions (namely zero field in the interior of the star, and zero net charge on 
its surface). Since the electric field is zero inside the star, the potential V s + V should be a constant on the surface of 
the star, hence Vf(i? s ) + V t (R„) = for £ ^ 0. For reasons already stated in section HV Al Vg(R) = Cgfe(R) outside 
the surface of the star. Using (fl~8|) for Vf , we can write 

eg e (R s )f e (b) + C e f e {R s )=0: 1^0. (42) 

This equation can be solved to obtain Ci for all £ ^ 0. 

Co can be obtained by imposing the condition that the net charge on the star vanishes, that is <7° = 0. To impose 
this condition, we should first apply Gauss's law to obtain an expression for the induced surface charges 5(0). Since 
the electric field inside the star vanishes, applying Gauss's Law simply gives (see appendix rA"|) : 

. f -R 2 S [e fe(b)g/(R s ) + C,f/(R S )] ^0 for 

47TCT = < 

( - R 2 S [eM/R 2 s b ~ C /R 2 ] for I = 0. 

In particular, using the fact that ct° = 0, we find that Cq = eM/b. Note that Vo(r) = 2e/b for r < r s . 



R 



dV(R,i 



R=R S 



dR 



(43) 



R=R. 



1. Evaluating the self force 



Now that we know all the coefficients CV, we have completely determined the regular field that would produce the 
self force. As in the previous section (insulating star), (|31)|) gives the self force on the charge, with the only change 



10 



- 




Insulating 


■ 


- 











/ conducting 
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FIG. 2: The self-force on the charge is plotted as a function of its position for the two cases considered in the text: when the 
star is an insulator (upper curve) and when it is a conductor (lower curve). 



being that the Ci are now given by ([42]) rather than (f33|) . So, the self force is 



F sc \f = Kelt = — fo ( b ) + e 2_^\ j , n ^ ) ft V>). 

e e=i 



MR.) 



(45) 



Observe that the result obtained for the self-force in (|45|) docs not depend on the interior metric of the star. 

To compare the results with that of the insulating star, we assume that the size of the star R s is the same as that 
of an extreme density star, R s = 9M/4. Once again, it turns out that each F^ elf can be expressed as 



M 
~~b 



elf b 2 

We again define x = M/b for convenience. The form of ( l (x) is given for I = 0, 1, 2. 



(46) 



(47) 



C 2 (*) 



108 



a; 4 (72 In 3- 160) \2 



In [1 - 2x] + x + x' 





'1 


"(- 




900 




(1080 In 3- 


1376) ^ 



In [1 - 2x] - x + x' 



3 9 2 
---x + 3x 



9 

r 



In [1 - 2x] - 6x 2 + x 3 



3x 



(48) 



In [1 - 2x] + x 3 + 3x 2 - 6x) 



(49) 



Whereas the r] e (x) are positive for all values of £, the C^( a; ) are negative for all I except I = 0. For large £, 
(fix) — » 0, for all x. The convergence of C, l (x) is also plotted in Fig [1] for various values of x. We again note that 
as x — > (b — > oo), we have F se if — > i 7 ^; a charge placed far away from this star, would feel the same force as 
it would feel if the star were replaced by an equally massive black hole. However, something interesting happens 
nearby. There exists a point, xq, where the charge feels zero self-force. This is an electrostatic equilibrium point. It 
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happens when YleLi C £ ( x o) = ~C°- This condition can be numerically solved, giving xq = 0.34 (b = 2.95M). When 
x < x (b > 2.95M), the charge would experience a repulsive force and when x > x (b < 2.95M), the charge would 
experience an attractive force. 

As an example, we shall numerically compute the force experienced by the charge placed at x = 0.42. 



-Fsclf — 72 



M 

T 



S,C £ (0.42) 



-19.96-^ 



M 
~~b 



(50) 



The closer we take the charge to the surface of the star, the more the attractive force increases. In the limit x — > 4/9, 
-Fseif — * — oo. This is expected because even in flat space, when a charge is placed close to a conductor, it experiences 
a huge attractive force. We should interpret the attractive force which dominates when the charge is close to the 
conducting star as a consequence of electrical polarization of the star and the repulsive force which dominates when 
the charge is far away from the star as due to interaction of the fields with space-time curvature. Fig [5] presents the 
self force experienced by the charge F se if (normalized with respect to F^) as a function of the position of the charge 
placed near a conducting and an insulating star. 

A simple consistency check can be performed in the limiting case when M — > 0, by comparing it to the well known 
flat space limit. In flat space, the self force on the charge e at a distance b from the center of a conducting sphere of 
radius R s and zero net charge can be calculated to be 



rifiat 
self 



b 2 



(51) 



In the limit M -> 0, (|19I20[> reduce to f e (R) -> 1/R e+1 and g e {R) -> R l . As a result (gSJ) reduces to 



e=i 



which exactly matches the flat space result ([STj). 



APPENDIX A: GAUSS'S LAW 



In this appendix, we apply Gauss's law to evaluate the charge density on the star's surface. We start by reviewing 
Gauss's law. Consider a spacelike 3-hypersurface E, bounded by a 2-surface 5E. Then Gauss (Strokes) theorem can 
be written as [l|| 

/ VuF^dLp = \l F^dS^. (Al) 

From Maxwell's equations, the LHS corresponds to — 4ir times Q, the total charge enclosed within the bounded region 
E. The surface element dS^ = —2n\ ll r ) v \- l fad6d(f), where n M is the future directed unit normal to the space like 
hypersurface E and 7„ is the unit normal to the boundary <9E. The vector 7^ should always be directed away from 
the surface <9E. The term y/adOdcf), represents the 2 dimensional surface area element. 

In this paper, since we are interested in the surface charge density of the star, we apply Gauss's law in a region 
very close the surface of the star. We choose the hypersurface E to be a constant time T (or t) hypersurface spatially 
bounded by the Gaussian surface <9E, which is made up of two small pieces of constant- R surfaces of coordinate area 
dOdc/), one just outside the star surface R = R s ([T]) and the other just inside the star surface r — r s ([8]). 

For the piece of <9E outside the star's surface, in the Schwarzschild coordinate system (T, R, 6, (j>) as in metric 
(HJ), n p = - 2M/R, 0,0,0) and 7^ = (0,1/^1 - 2M/R,0,0). The only nonvanishing component of dS^ is 

dSrn = R 2 s'mdd8d(f>. For the piece of 9E inside the surface of the star, dS^ v is obtained from the unit normals to 
constant-t and constant-r surfaces with respect to the star's interior metric. With respect to the star's coordinate 
system (t,r,6,4>) as in metric (|5J|, = {—(f3 2 + ar 2 ), 0,0,0) and 7^ = (0,-1,0,0). Here, the only nonvanishing 
component of dS^ u is dSt r = — (/3 2 + ar 2 )r 2 sin 8 d9d<fr. 

By applying Gauss's Law, (|A1|) . on E, we obtain, 

- 47^(0, <{>) sin 9dBd4 = g TT 9 RR F TR dS TR \l u X C + 9* Q rr F tr dS tr \™* e . (A2) 
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Here, a(6, 4>) is the surface charge density and a(8, <j>) sin ddddtfi is the total charge contained within that small Gaussian 

surface. The axial symmetry of the problem considered in this paper ensures that cr((9, 0) is actually just cr(6). We 
have adopted (|A2j) in section llVl of this paper. For the case of conducting star (section [iVB[l . we note that the field 
inside the star vanishes and hence the second term in the RHS of (|A2|) vanishes. 



1. Conformal invariance of Maxwell's equations 



As already indicated, we can ignore the conformal factor in (j4]) and use the metric ([8]) when evaluating the integral 
i J F^dS^v on the surface just inside the star. This is justified as long as the source {J^\J—g) is conformally 
invariant, in which case: 

• the LHS of (|Aip (that is, the total charge contained within a coordinate volume) is invariant under a conformal 
transformation, 

• Maxwell's equations have the same solution for the vector potential A^, and hence the same solution for F^ Vl 
in any conformally related metric. 

Note that, under a conformal transformation g^ v — ► VP'g^ (g^ — ► Q, g^ u ), so we also have — > On^, 7 M — * ^7^, 
•v/ — g — s- fl 4 ^/ — 5 and y/a — > ft 2 \fa . Hence, 

dS^ n 4 dS^. (A3) 

Also note that, if the source (J^yf—g) is conformally invariant, Maxwell's equations (|lfjp require F M „ to be conformally 
invariant. With contravariant indices F^ v transforms under a conformal transformation as 

aft Q pv 1 

* -g g ^2 n2 *a0 ^ ■ lA4j 

Cleary, (|A3[) and (|A4[) imply that (|Aip is invariant under a conformal transformation. 
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